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Correlation Functions of Harish-Chandra Integrals 
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A. Prats Ferrer^^, B. Eynard* 0, P. Di Francesco^^, J.-B. Zuber^^ 



Abstract 

The Harish-Chandra correlation functions, i.e. integrals over compact groups of invariant 
monomials Utr [XP^nV^Q^XP'' ■ ■ ■) with the weight exp tr (XfiFfi''") are computed for the 
orthogonal and symplectic groups. We proceed in two steps. First, the integral over the compact 
group is recast into a Gaussian integral over strictly upper triangular complex matrices (with 
some additional symmetries), supplemented by a summation over the Weyl group. This result 
follows from the study of loop equations in an associated two-matrix integral and may be viewed 
as the adequate version of Duistermaat-Heckman's theorem for our correlation function integrals. 
Secondly, the Gaussian integration over triangular matrices is carried out and leads to compact 
determinantal expressions. 
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1 Introduction 



In the study of matrix integrals [B [21 [3l H] , one frequently encounters integrals of the form 

Z^= [ rffie-*^(™") (1-1) 



over some compact matrix group G, with X and Y two given matrices. By the left and right 
invariance of the Haar measure dfl, this integral is invariant under 

X iliXn^^ , Y fisl'^^^^ , (1-2) 

and is thus insensitive to the choice of the representative of the orbits of X and of Y under the 
(adjoint) action of the group. This enables one to bring the matrices X and Y to some canonical 
form, as we shall see below. 

The case of reference is the so-called Harish-Chandra-Itzykson-Zuber (HCIZ) integral [HIE], 
where the integration is performed over the unitary group Q G U(n) and X and Y are two 
(anti)Hermitian matrices. By the previous argument, we may with no loss of generality assume 
that X and Y are two diagonal (anti)Hermitian matrices of size n, X = diag(Xj)j=i^... and 
likewise for Y. 

where 

A{X) = ll{X,-X,) (1-4) 

i<j 

is the Vandermonde determinant of the eigenvalues Xi of X and likewise for A{Y); Y^ = 
diag(y7r(j))j=i,... Further examples are provided by more general Harish-Chandra-type inte- 
grals, where X and Y live in (a matrix representation of) the Lie algebra g of G [5j. For example, 
G = 0{n), X and Y anfeymmetric real matrices of size n. In all these cases, explicit formulae 
are known, following from a diversity of methods, see below and for example [7] for a review and 
references. 

It is desirable to extend these formulae to the "correlation functions" of the integral f ll-ip . i.e. 
integrals of the form 

J dnF{x,nYn~^)e~''i^^^^'') (1-5) 
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with F invariant under (ll-2p . Such correlation functions provide a deeper probe of these integrals 
and, in a physical context, give often access to quantities of interest. They also act as generating 
functions of integrals of the form 

j dn ■ ■ ■ ■ ■ ■ e"*' ^'"'"''''"^ (1-6) 

i.e. of moments of Q and fi^ with the Harish-Chandra weight. 

Kogan a/ [8], Morozov [9], and Shatashvili [10] made some attempts at computing correlation 
functions of the HCIZ integral, i.e. for the unitary group. Morozov's formula may be recast into 
a very compact expression [TT] but it is only good at computing correlators quadratic in fi, whilst 
Shatashvili's formula allows in principle to compute all correlators, but is not of easy use. In paper 
|12] two of us have shown how to recast the computation of correlators of the type fll-5p for the 
unitary group into a totally different setting. The method consists in two steps. In step one, the 
integral is rewritten as a sum of integrals over upper triangular complex matrices. The formulae 
of [12] are in some sense a generalization of Morozov's, and allow to compute all correlators for 
the U{n) group in a very simple formula. The initial observation was that a Gaussian integral 
and its polynomial moments in an hyperplane of dimension d of M„(C)^ does not depend on 
the hyperplane (up to multiplication by a constant jacobian), and thus, one can use either the 
hyperplane Mi = M|, M2 = or the hyperplane Mi = M^, which have the same dimension. In 
the first hyperplane, after diagonalisation of Mi and M2, the integral separates into a radial part 
and an angular part proportional to moments of U(n) with the HCIZ measure, while in the second 
hyperplane, after Schur decomposition of Mi, the integral separates into a radial part (identical 
to the one in the first hyperplane), a trivial angular part, and a Gaussian integral over triangular 
matrices. As a result, the authors of [12] were able to identify all moments of the HCIZ integral 
with a Gaussian integral over complex strictly upper triangular matrices. The formula reads: 

/ df/F(X,[/Ff/t)e-*^(^^^^') = 
^(_l)-e-*'-™ f dTF(X + T,i; + Tt)e-*'-(^^') , (1-7) 



A(X)A(F) 

for X and Y two real diagonal matrices. 

In a second step, the Gaussian triangular integrals in the right hand side were computed in 
[12], using Wick's theorem. The computation can be performed explicitely due to the nilpotent 
properties of T, which ensures that most Wick's pairings actualy vanish. The computation is most 



easily done by recursion on the size n of the matrix, i.e. by integrating out the last column of T. 
An appropriate basis of all possible polynomial moments F was introduced in [12], and in that 
basis, it was found that: 

» n 

/ dTF{X + T,Y^ + T^)e-''i^^') =llMiX,,Y^^,)) , (1-8) 

where F and Ai{x,y) are matrices of some size R\. The universal matrices Ai{x,y) have many 
remarkable properties, in particular they commute with one another 

[M{x,y),M{x',y')] = . (1-9) 

The purpose of the present article is to generalize the computations of [12] to other classical Lie 
groups. Specifically, we address the computation of f ll-5p for G the orthogonal group 0(n) or the 
symplectic group Sp(2n), with X and Y in the Lie algebra of those groups. First, we relate H-C 
correlators over those groups to Gaussian integrals over some set of triangular matrices, then we 
compute the latter Gaussian triangular integrals using an appropriate basis, and finally we find 
that the result can again be written as products of the same matrices Ai{x, y) which appeared for 
U(n). Our main results are stated in Theorems 14. 3[ [573] and [6A] below and in section [3, (eq. fl7-2p . 
eq. ( 17-3P and eq. (17-40 ). The results for the unitary, orthogonal and symplectic groups may be 
expressed in a unified way, in terms of the Weyl group W, Borel subalgebra b and positive roots 
a, in the following 

Theorem 1.1 

/ dnF{x^,nY^n-') 6-*^^"^^"^^-^) 

Jg 

. f,+ tr {Xw(Y)) r . , 

= TT / dTF{^X + T,^wiY)+T^)e-'^i^^') (1-10) 

for any polynomial function F and for some F-independent constant c; here X*^ and Y"" are taken 
in a Cartan subalgebra, and should thus be thought of as aniz-Hermitian matrices with extra 
symmetries depending on G, (see sect. 2 for more details), while iX and iY are the purely imaginary 
diagonal matrices with the same eigenvalues as X"- and Y"-. In that form, the derived Borel 
subalgebra n+ is made of complex strictly upper triangular matrices, also subject to symmetries. 
It is thus natural to expect these results to extend to any simple compact group G, see Conjecture 
8.1. 
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It is our hope that these results should provide a new insight on the common features of all 
these integrals. 

Our paper is organized as follows. In section [2l we review the known results by Harish- 
Chandra and Duistermaat-Heckman and set up the notations. In section [3], we show how Gaussian 
integrals over two matrices with reality properties, either antisymmetric real or antiselfdual real 
quaternionic, may be equated to Gaussian integrals over one complex matrix constrained by some 
symmetry requirements. This is established by use of loop equations, on which we provide details in 
AppendixlHl Section [Hand [5] then show how separation of the angular variables by diagonalization 
or Schur decomposition leads us to the desired integrals, which are thus related to integrals over 
complex triangular matrices (with additional symmetry requirements). Section [6l supplemented 
by Appendices O and [El is devoted to the actual computation of these integrals over triangular 
matrices, by means of a recursive method using a diagrammatic method. The final expressions are 
displayed in section [TJ while section [8] contains our concluding remarks and suggestions of further 
directions worth exploring. Two other appendices make our notations explicit on quaternions 
(App. or give additional details on the calculation of some Jacobians (App. [Cj. 

2 Overview of known results 
2.1 The Harish- Chandra theorem |[5] 

Following Harish-Chandra, for G a compact connected Lie group, we denote by Ad the adjoint 
action of G on its Lie algebra g, by (X, Y) the nondegenerate invariant inner product on g, which 
we take to be the trace of the product XY in our matrix representation, by f) C the Cartan 
subalgebra, and by a{X) the linear action of a root a on X G f). If X and F G () 

A(X)A(r) / exp-(X,Ad(l])r) = e(^)exp-(X,^(F)) , (2-1) 

where w is summed over the Weyl group W, e(w) = (— l)'*'^"'^ ^{w) is the number of reflections 
generating w, and 

A(X) = n«W, (2-2) 

a product over the positive roots of g. A(X) may be called a generalized Vandermonde de- 
terminant, since in the case of f/(n), it reduces to eq. fll-4p (and (12-11) reduces to (I1-3P ). while 
the expressions for the orthogonal and symplectic groups will be given below. The constant 
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<7r, 7r> in the right hand side of (12- ip is computed as follows. Write all the positive roots in 
an orthonormal basis Ei of root space, i = 1, ■ ■ with i the rank of g. Regard vr = no>o'^ 
as a polynomial in the positive roots and expand it on symmetrized tensor products of the e, 

^ = Em,>oP("^i'"^2,--- ,"^£)£r ■■■^r- Then <7r,7r>= Em,>o(p("^i' "^2, ■ ■ ■ , IlLi "^i- 

For = su{n), one finds <7r, 7r>= nj=i j'' while the expression for the other classical groups will 
be given below. 



2.2 The Duistermaat-Heckman theorem 1 1131 



The Duistermaat-Heckman theorem states that if is a symplectic manifold, invariant under 
a U(l) flow generated by a Hamiltonian H, then for the integral Jj^ e*^*, the stationary phase 
approximation is exact: the sum of the values of the integrand at its critical points, weighted by 
the Gaussian ('one-loop') fluctuations around them, gives the exact integral. 
For X and Y in the Lie algebra, consider the integral eq. (11-11) 

Z = [ rffie-*^(™") . 
Jg 

In such integrals, we first pick a convenient representative of the orbits of elements of the Lie 
algebra under the adjoint action of G. A theorem of Cartan asserts that any element of the 
Lie algebra is the conjugate (under the adjoint action) of an element of the Cartan algebra |15j . 
Thanks to this theorem and to the left and right invariance of the Haar measure dQ, one may 
always assume that X and Y lie in the Cartan subalgebra i). This assumption matches that 
of Harish-Chandra's theorem. Moreover the integration is then reduced to G/T, T a maximal 
abelian subgroup (Cartan torus) commuting with Y, or alternatively, the integration is carried 
out on the orbit of Y under the action of this quotient. This is a symplectic manifold, to which 
Duistermaat-Heckman's theorem applies |14j . 

We thus first look for the critical points of the 'action' tr {XQYQ~^) when Q G G/T. In other 
words, we look for solutions in Q of 

5tr {xnYn~^) = tr {6m~^[x,nYn~^]) = o , 

Since A := is arbitrary in 0\f), this implies that the component of [X, in 0\f) 

vanishes. On the other hand, the component of [X, fiFfi"^] in () also vanishes, since if B := 
[X,nYn~^] were in [), then tr {Bf = tr {B [X,fiFfi-i]) = tr {[B,X]nYn~^) = since X and 
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-B G [) commute. We thus conclude that B = 0, i.e. that 

[x,nYn~^] = o . (2-3) 

The critical points are thus the points Qc G G/T such that (12-31) is satisfied, which for generic 
X G f) means QcY^c^ ^ ^-S- takes the element F G P) to an element QcY^c^ ^ ^- 
denote by W the normalizer of the Cartan torus T quotiented by T, the previous discussion has 
just proved that the critical points Qc of the action are in one-to-one correspondence with elements 
w of the group W. The group W is known to be the Weyl group of G ([I5], Prop. 15.8). In the 
sequel, we denote = fi^Ffi"^ for w G W. 

At this stage, Duistermaat-Heckman's theorem thus tells us that 

/d(]e-*^(™"') = V / ciAe-[*^(^'^"^"^"")]^ (2-4) 

where [■ ■ ■]2 means that we retain only up to the quadratic terms in the expansion in powers of 
Aeg\i). 

The final step in the application of Duistermaat-Heckman theorem is thus to compute the 
second order variation of the action at one of these critical points. For Q = e^, A G g\i), 

- tr [Xe^Y'^e-^] = - tr {XY"") + ^ tr {[A, X] [A, F""]) + o{A^) . (2-5) 

We then have to carry out the Gaussian integration 

"■^.^^|tr {[A,X][A,Y-]) 



over the (i-dimensional vector A. This (real) dimension d = dim g — dim P) = 2r is even and equal 
to the number of roots of G. We now expand A, X and Y^ in the standard basis A = AaEa, 
X = J2i ^iHi and Y^ = J2i ^^"Hi and use the standard commutation relations and traces 
tr {HiHj) = 6ij, tr (EaEp) = 6a+p,o to get 

tr {[A,X][A,Y-]) = J2 A^ApX,Y;"a^'^(3^^hT (E^Ep) = - ^ l^'^a^^') 

a,l3,i,j i j 

i.e. 

tr ([A,X][A,F-]) = -^A,A_,«(X)«(F-) 

a 

with a sum over positive and negative roots. This quadratic form has a signature (-!-'', — ^), and 
upon a suitable contour rotation, the integration over A yields 

/ d'^A i[AX][AY-]) ^ ^Q^stant 

J n.>o«(^)«(>^"') ■ 

Putting everything together, we see that we have reconstructed the Harish-Chandra formula. 
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2.3 Explicit formulae 

It is of course a good exercise to repeat these steps and to write explicit expressions for each 
of the classical groups U(n), 0(n) and Sp(n). The result for U(n) is well known and has been 
recalled above. We shall content ourselves in giving the final result for the two latter cases. In 
the orthogonal case 0(?7,), we have to distinguish the n = 2m and n = 2m + 1 cases. In the even 
case, G = 0(2m), we take the X and Y matrices in the block diagonal form 

Then the critical points Qc are the product of a permutation r of the m blocks of F by a diagonal 
matrix of signs diag(ti, . . . , tm), where tj = ±16.2, or in other words, the set W is x Z^. Note 
that W is larger than the ordinary Weyl group of Dm = so(2m) type, which is 1^ = &m x Z™~^: 
this is because changing the sign of one Yj, say Yi, is performed by conjugation by a matrix made 
of 2 X 2 blocks, Q = diag(cri, Id2, ■ ■ ■ ,Id2), which is in 0(2m) but not in S0(2m). As a result, 
only an even number of signs may be changed in the latter case, whence the factor Z™^^ in the 
Weyl group. For the 0(2m) group that we consider here, we thus have 

Zioi2rn)) ^ y -r^-;ry^ = const, y e ^^'^\, , ^ 

w={t,{U}) ^ ^ ^ > tG6^ \ } \ } 

rnn^t det(2 cosh (2X,F,)),,,=i,.,^ 
" A(X)A(F) ' 

where Sr is the signature of the permutation r, and 

A(X) = n(^- - • (2-8) 
For n = 2m + 1, the calculation proceeds along the same line. We write 

'0..,„,j») o),,.,„,„-°) ■ 

The critical points Qc are again the product of a permutation r of the m blocks of 5 by a matrix 
of signs, tj = ±Id2, j = 1, ■ ■ ■ , m, and 

^(0(2^+1)) ^ y ; ' I = const, y e^-^^^^——^^ L 

Mu}) ^(^)^(^'") A(X)A(r) 

_ det (2 sinh(2X.y, )).,=,,.., ^ 

- const. A(X)A(r) • 
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with now 

m 

A(X) = n(^- -^DII^- (2-11) 

Finally for the symplectic group Sp(2m), it is convenient to use quaternionic notations for 
matrices, i.e. to regard the matrix elements as quaternionj^l, with coordinates in the standard 
quaternion basis, eo = l; ej = —1 , i = 1, . . . ,3 , 6162 = 63 ; alternatively, the matrices may 
be regarded as made of 2 x 2 blocks written in terms of the identity matrix Id2 and of Pauli matrices 
a (with the identification Cq ^ Id2 , ej ^ ""^crj) j = l!2,3). The Lie algebra Cm of Sp (2m) is 
thus generated by quaternionic real and antihermitean (also called antiselfdual quaternionic real, 
see appendix m x m matrices X, Xij = X^- + Xij e, X^- G M, X = —X'^. Consider the Cart 



an 



algebra generated by the m matrices diag(Xje2), j = 1, ■ ■ -m. We thus take our matrices X and 
Y of that form 

X = diag{Xje2)j=i,..,m , Y = diag{Yje2)j=i,.. ,m ■ (2-12) 

Then, the critical values Qc are again the product of a permutation r of the m blocks of F by a 
diagonal matrix of signs, tj = ±1. This leads to 

, = const. V e ^^^ ^ ^ (2-13) 

_ deti2sinhi2X,Y,)),,=,,...,m . . 

- const. A(X)A(F) ' 

with the same expression for A(X) as in fl2-lip . Thus Z'^^p^^™^^ has the same form as the integral 
over 0(2m + 1). 



'We refer the reader to appendix [A] for more details on our notations on quaternions. 
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2.4 List of notations 

For the sake of the reader, we hst hereafter the non standard notations in the order they appear 
in the text. 

A(X) Vandermonde determinant and generahzations (11-41) . fl2-2p . fl2-8p . (12-1 ip 

An n X n real antisymmetric matrices (13-11) 

J antidiagonal identity matrix (13-31) 

J An n X n J — antisymmetric complex matrices sect. 13.1.21 and (l3-4p 

QAm m X m real quaternionic antiselfdual matrices sect. 13. 2. H and ( ]3-7p 

J antidiagonal symplectic matrix (13-111) 

JA2m 2m X 2m J — antisymmetric complex matrices sect. 13.2.21 and (13-131) 

-D^(]R) real 2x2 block — diagonal antisymmetric n x n matrices sect. 14.11 

M„(C) n X n complex matrices sect. 14.21 

Tn n X n strictly upper triangular complex matrices sect. 14.21 

D„(C) n X n complex diagonal matrices sect. 14.21 

U'^(?7,) twisted orthogonal matrices (14-51) 

n X n strictly upper triangular fTHI) 
J— antisymmetric complex matrices 

D'l^iC) , -D^(]R) n X n complex, resp. real, J— antisymmetric diagonal matrices (14-71) , (I4-18P 

jjaRr^ m X m real quaternionic diagonal matrices ^^^^ ^-^^ 

™ with elements proportional to 62 

U'^(2m) twisted symplectic matrices (15-61) 

2m X 2m strictly upper triangular ((5^ 
^"^ J— antisymmetric complex matrices 

3 Analytical continuation for two-matrix integrals. 



In this section, we follow the same strategy as used in [12] for the unitary group: the integrals 
of interest ( 11-5P are regarded as the "angular part" of two-matrix integrals over the classical Lie 
algebras so(2m-|- 1), so(2m) and sp(2m), and the latter may be analytically continued to integrals 
over complex matrices with special symmetries. 
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3.1 Real antisymmetric two matrix integral and complex J-antisymmet 
matrix integral. 

3.1.1 Real antisymmetric two matrix integral. 

Consider first the set An of nxn real antisymmetric matrices and consider the measure on An x An 

d/i(Ai,A2) = e-*'-('^^?+'?^'+^^^^^)dAidA2 

dAk = n ; k = l,2 (3-1) 

l<?<j<n 

Then the real antisymmetric two matrix partition function and the associated correlation functions 
are defined as 

Z2RA = / dfi{Ai,A2) 

{F{A„A2)),^^ = ^ [ F(Ai,A2)d/i(Ai,A2) 

using the measure dfi{Ai, A2) given in eq. fl3-ip . The partition function is the product of n(n — 1)/2 
uncoupled and equal integrals over the pairs of matrix elements {{Ai)ij, (742).^), i < j. Each 
integral, of the form f dxdy exp{{x,y)Q{x,y)^}, Q = i'^^ ^ ) , is absolutely convergent if the 
real part of the quadratic form {x,y)Q{x,y)^ := aix^ + a2y'^ + 2jxy is negative definite, which 
holds true if ReaiRea2 — (Re 7)^ > and Reai, Re 0^2 < when x and y are integrated over the 
real line. Then the partition function is easily computed to be 



Ti { n — 1 ) 



Z.,. - (Tj) ' (3-2) 

where 6 = a\a2 — 7^. Likewise, for polynomial Fi^A^^A^)^ the correlation function {F)^^^^ is by 
Wick theorem a polynomial in the matrix elements of the propagator Q"^, namely ^ , ^ and ^. 

3.1.2 Complex J-antisymmetric matrix integral 

Define now the nxn antidiagonal matrix J = J^^ 

/ ■■■ 1 \ 



J 



(3-3) 



\1 ■■■ Oj 
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Any matrix M with the property JM^ = —MJ is said to be J -antisymmetric. Such a matrix is 
antisymmetric with respect to the second diagonal, i.e. Mij = —Mn+i-j^n+i-i 

( Mi,i Mi,2 ■■• Mi,„_i \ 
M2,i ••■ ••■ -Mi„_i 



V -M„_i,i 



-Mi,2 

-M2,i -Afi,i y 



and in particular, Mij = whenever i + j = n + 1. 

On the set of complex J-antisymmetric matrices, we consider the measure 



d/x(M) 
dM 



. tr (i^M2 + ^Mt2+7MtM) 



dM 



JJ dReMi J dImM, 

i+j<n+l 



(3-4) 



Then the complex J-antisymmetric matrix partition function and the associated correlation 
functions are defined as 



n.jA 



d/i(M) 



IJA 



dfi{M)F{M, M^) 



ZlJA JJA„ 

using the measure d/i(M) given in eq. fl3-4l) . The partition function is again the product of the 
n{n — l)/2 uncoupled and equal integrals over the complex independent matrix elements Mij, 
i + j < n + 1. It is absolutely convergent if Reai, Rea2 > 0, Re7 < and Re7^ > Reaia2 and is 
then given by 



ZijA 



TT 



2v^ 



n(n—l) 



(3-5) 



with S = axa2 — 7^ as before. For polynomial F(M, M^), the correlation functions {F)^jj^ are 
again given by polynomials in ^ , ^ and ^. 

3.1.3 Analytic continuation. 

The two families of integrals just studied have close connections, even though their original domains 
of convergence may not overlap. The first trivial observation is that in both cases we have the same 
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number of variables, as already manifest in the computations of the partition functions, namely 
n{n — 1) integration variables in both integrals. The second and more important observation for 
our purpose is that both integrals share the same loop equations. This will be proved in appendix 
iBl Third, as already stressed above, the correlations of polynomial invariant functions both in 
An X An and in J An are polynomials in the variables ^ and j, and thus analytic functions. 
From all these observations, we formulate the following 

Theorem 3.1 The polynomial correlation functions of the two real antisymmetric matrix integral 
are equal to the correlation functions of the complex J -antisymmetric matrix integral in the sense 
of analytic continuation, i.e. 

{F{Au A^)),^^ = {F{M, Mt))^^^ . (3-6) 

Proof: Note that the loop equations given in appendix [B] are in fact recursion relations on the 
polynomial degree of the correlation functions. Then the fact that polynomial invariant correlation 
functions are polynomials in ^ and j and the fact that the loop equations and their initial 
condition (namely (1) = 1) are the same for both integrals imply that the polynomials generated 
from the recursion are the same. □ 

Although the correlation functions are not originally defined in the same region in parameter 
space, the fact that they are polynomials allows one to analyticly continue them and to identify 
them. 

3.2 Real quaternionic antiselfdual two matrix integrals and complex 
J-antisymmetric matrix integral. 

In this section we consider another pair of matrix integrals, related to the symplectic group, for 
which similar considerations hold true. 

3.2.1 Real quaternionic antiselfdual two matrix integrals 

Consider first the set QAm of real quaternionic antiselfdual (antihermitian) mxm matrices, whose 
definition has been recalled in sect. 12.31 and appendix lAl 
On QAm X QAm, we consider the measure given by 



vi<j a=0 / \j=l a=l 
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where tro (• • • ) = 2Retr (...) is a scalar (while tr (...) is in general a quaternion number, see 
appendix Rl) . The quadratic form in this 'Gaussian' measure is thus 

3 

-tro (yg? + yQ^ + 7QiQ2) =2 J2 E ("1(^1 )1 + "2(^2^- + 27(Q?).,(Q?).,) 

l<i<j<rn a=0 

m 3 

+ E E + «2(Q2)?. + '^liQiUQ^h) ■ (3-8) 

i=l a=l 

The real quaternionic antiselfdual two matrix partition function and the associated correlation 
functions are defined as 



'2QA 



I MQ1.Q2) (3-9) 



^20 A JOAr, XOAr, 



^2QA JQAnxQA 

The partition function is readily computed to be 

and once again correlation functions of polynomials in Qi,Q2 are polynomials in and |. 

3.2.2 Complex J-antisymmetric matrix integral 

We now introduce a 2m x 2m matrix J = —J~^ of the form 

written in terms of J defined above in f l3-3p . Any matrix M with the property JM'^ = —MJ 
is said to be J-antisymmetric. Such a matrix possesses a peculiar symmetry with respect to the 
second diagonal: we can write it as 

f A B 
^=(c D 

where A, 5, C, D are m x m matrices satisfying 

A = -JD^J ■ JB^ = BJ ; JC^ = CJ (3-12) 

Thus, under the reflection with respect to the second diagonal, B and C are invariant, while A 
and —D are exchanged. 
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On the set JA2m of complex J-antisymmetric matrices we consider the measure 

d/i(M) = e-^=^(^*^^+^^^^+^^^^'OdM 

dM = Yl dReMij dlmMij . (3-13) 

i+j<2m+l 

Then the complex J-antisymmetric matrix partition function and the associated correlation func- 
tions are defined as 

ZijA = i MM) (3-14) 

'^IJA JjA2m 

using the measure d/i(M) given in eq. fl3-13p . The J-antisymmetric partition function reads 

and once again, correlation functions of polynomials in M and M ^ are polynomials in the param- 
eters f , f , and f 

3.2.3 Analytic continuation 

Again the observations made in sect. 13.1.31 extend to this case. The two matrix integrals of 13.2. H 
13.2.21 have the same number of integration variables equal to 2m(2m + 1), they satisfy the same 
loop equations (see appendix [B]), and their correlation functions of invariant polynomials have 
polynomial dependence on ^ and |. 

These observations allow us to formulate an analogous analytic continuation theorem for these 
two matrix integrals 

Theorem 3.2 The polynomial correlation functions of the two real quaternionic antiselfdual ma- 
trix integral are equal to the correlation functions of the complex J-antisymmetric matrix integral 
in the sense of analytic continuation, i.e. 

{FiA,, A2)),Q^ = (F(M, M^))^j^ . (3-16) 

Proof: The proof goes exactly as the one in theorem 13. 1[ □ 
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4 Correlation functions over the orthogonal group 

In this section we exploit the relation found in theorem 13.11 by performing a separation between 
"angular" and "radial" variables of matrices in the two sides of equation eq. (13-61) . 

4.1 Block-diagonalization of antisymmetric matrices 

We first consider the case of antisymmetric matrices and of the orthogonal group 0(n) equiped 
with its Haar measure dO (normalized to / dO = 1). 

As recalled in sect.l, Cartan's theorem asserts that any antisymmetric matrix A may be 
brought to the block diagonal form (l2-6p or (l2-9p by an orthogonal transformation of 0(n) (for 
this standard result, see also [161 [HI [IS]). Denote by -D^(M) the set of such real block-diagonal 
antisymmetric n x n matrices, with the Lebesgue measure: 

m 

dX:=l[dX,. (4-1) 

i 

By an abuse of language, we shall refer to the Xi as the "eigenvalues" of A. 
In the new variables {0,X}, the Lebesgue measure in An reads 

dA = Jac^ A\X) dO dX (4-2) 

where the Jacobian is (see appendix [O for details) 

Jac«= ^^^Wf^ (4-3) 
" I ^"^ 2- ifn = 2m + l 

We recall that A(X) takes two different forms (12-81) and (12-111) depending on the parity of n. 

This decomposition is unique up to a permutation of the m "eigenvalues" , a change of signs of 
each eigenvalue independently, and a multiplication of O by a 2 x 2 block-diagonal matrix whose 
diagonal blocks belong to O (2). In other words, A = 0X0^ establishes a mapping between An 
and O (n) x Z}"(R)/ (O (2)™ x (5^ x Z^). This overcounting has already been taken into account 
in Jac^. 

4.2 Schur decomposition of complex J-antisymmetric matrices 

A less standard result, (see flG l fT7 l fT8l fT9] for instance), is that any complex matrix M G M„(C) 
can be written as: 

M' = U'{Z' + T')U'^ (4-4) 
15 



where f/' G U (n) is a unitary matrix, T' E Tn a strictly upper triangular complex matrix and Z' G 
D„(C) a complex diagonal matrix. We can apply this Schur decomposition to a J-antisymmetric 
matrix. This will induce further constraints on the unitary and triangular matrices. 
Define U"^(?7.) to be the subgroup of U (n) satisfying the condition 

[/-I = f/t = JU^J (4-5) 

with the induced normalized Haar measure. We will call these matrices twisted orthogonal matrices. 
Define also to be the set of n x n strictly upper triangular J-antisymmetric complex matrices, 
with the Lebesgue measure: 

dr:= Y[ dReTijdlmTij (4-6) 

i<j 
i+j<N+l 

and D'l{C) to be the set of n x n complex J-antisymmetric diagonal matrices 

Z = diag{Zi,- ■ ■ ,Zm,-Zra,- ■ ■ ,-Zi) or diag{Zi, ■ ■ ■ , Zm,0, -Z„^, ■ ■ ■ , -Zi) , (4-7) 
depending on the parity of n, with the Lebesgue measure: 

m 

dZ =Y[ dReZi dlmZi . (4-8) 

i 

Finally we define for these matrices of /^^{(C) 

AiZ)^Y\iZf-Z]).y ;[ ^llll^ . (4-9) 

i<j "^^^^ 

With these notations one can prove 
Proposition 4.1 Any J-antisymmetric matrix M may be written as: 

M = U{Z + T)U^ (4-10) 
where U G U'^(n), T E and Z G D:^{C). The Lebesgue measure in J An is then: 

dM = Jac^' \A{Z)\^ dU dT dZ (4-11) 

where the Jacobian is 
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Proof: Consider the Schur decomposition (14-41) of the matrix M. Noticing that 

det (A - M) = det (A - JMJ) = det (A + M) (4-13) 

we immediately see that the non- vanishing eigenvalues come in pairs (A, —A). By a possible 
redefinition of U, we may always order the eigenvalues in a J-antisymmetric diagonal form Z as 
in (14- 7p . The constraints on U and T follow from the J-antisymmetry of M and Z. The measure 
can be computed using the same method as in the appendices of [16]. □ 

This decomposition is unique up to a permutation of the m different eigenvalues, to changes 
of sign of the m eigenvalues and to multiplication of f/ by a diagonal matrix V G l]'^{n) whose 
elements are on the unit circle. In other words, M = U{Z + T)U'^ provides a 1-to-l mapping 
between JA(C) and \J\n) x T;[ x Di{C)/{\J (1)™ x 6„ x Z'2"). The overcounting is included in 

dHU). 

4.3 Orthogonal and triangular matrix integrals 
4.3.1 Radial and angular integrals 

Consider the block-diagonal decomposition of the real antisymmetric matrices and the Schur 
decomposition of the J-antisymmetric complex matrices. Using these decompositions we will 
rewrite both sides of eq. (13-61) . 

Theorem 4.1 A matrix integral over AnXAn can be decomposed into a "radial" and an "angular" 
part using the block-diagonal decomposition A = 0X0^ . 

[ di2{Ai,A2)F{Ai,A2) = (Jaa^)' / dXdyA^{X)A^{Y) 

J An X An iz)f^(R)xD°(R) 

^-tr{^xH^Y^) j^^^^ rfOF(X,OFO^)e-^*^(^'^^'^") (4-14) 
with the notations of ^2-8\) and ^2-ll\) . 

Proof: The theorem follows from the results of section 14. 1[ □ 

Notice that one of the two orthogonal matrices decouples and so this part of the integral gives 
1. The remaining orthogonal integral represents the relative angular variables. 
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Theorem 4.2 A matrix integral over J An can be decomposed into a "radial", an "angular" and 
a "triangular" part using the Schur decomposition M = U{Z + T)W . 



[ dfi{M)F{M,M^) = Jac'^' f (iZ|A(Z)| V*"" (^^'+'^^*') x 

J J An JdJ{C) 

[ dTF{Z + T,Z* + T^)e^^ (^^^) (4-15) 



-7tr (Z'Z) 



with the notations of ( [^-g| ). 

Proof: The theorem follows from the results of section I4.2[ □ 

Notice that the twisted orthogonal matrix decouples and so this part of the integral gives 1. 
Only the triangular and radial parts remain. Notice also that the measure for the triangular part 
factors out from that of the radial part and only a Gaussian measure remains for the triangular 
part. 

4.3.2 Relating integrals over orthogonal and triangular matrices 

In this subsection, we relate the HC integral over the orthogonal group O (n) 

jOW^^ f dOe-^*'•(™"°")F(X^OF'^0^) , (4-16) 

Jo{n) 

with X", Y°- G Z)^(]R), to an integral over complex upper triangular matrices of T^(C). Note 
first that X^*-"'' is a completely symmetric and even function of the "eigenvalues" Xi and of the 
Yi, i = 1, ■ ■ ■ ,m = \n/2\. This is because any permutation or sign changing matrix acting on 
either X"- or may be absorbed into a redefinition of the orthogonal matrix O. In contrast, the 
integral over complex triangular matrices will have to be symmetrized by hand. 

To obtain the desired relation between HC-type integrals over the orthogonal group and inte- 
grals over the triangular matrices of T^, we shall follow the same steps as in [12], in particular of 
Lemma A.l there, which asserts that for any polynomial u in two variables, one has the relation 



(4-17) 

where we have one complex variable integration on the left hand side and two real variables on 
the right hand side. This relation may be promoted into the following equality between integrals 
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over diagonal matrices 



(4-1^ 



We now apply theorems 14.11 and 14.21 to the two sides of (13-61) 



{F) 



2RA 



/ dAidA2e-*'-(^^'+'?^2+7^i^2)^(^^^^^) (4_19) 

Z2RA J An X An 



Z2RA Jdz 
X 



(JacJ' /■ dXMF"A2(X'^)A2(F'^)e-*^('^^"+^^") (4-20) 

>/ds(IR) xD«(M) 
Join) 

= i^F)uA = ^l dMe-*'-(^^-^'+'?*^''+^^^*^')F(M,Mt) (4-21) 

ZlJA JjAn 

= ^^'"'^^"(Q ^ / dZ|A(Z)|V*=^(^^^+^^"-^^^*) (4-22) 

x[ dT e-^''(^^')F{iZ + T,-iZ* +T^) . (4-23) 

In the last line, we have performed a change of variables Z ^ iZ for reasons that will appear 
soon. We then apply (14-181) to get 

{F)uA = ^^fz^m f dXdFA(X)A(F)e*^(^^^+^^^-^^^) 

X [ dTe-^^'^^^^)F{iX + T,-iY + T^) (4-24) 



^^'"^^^"^^^ / dXdFA2(X)A2(y)e*=^(^^^^-^^^) 

ZijaZd^ (r) X (R) 2'"m! J j ^ (r) 

E^TYYaT^/ ^Te--(-^)F(.X + T,-.r(.,,+Tt) (4-25) 



X 



In the last line, we have symmetrized the integral over triangular matrices for the reason 
explained at the beginning of this subsection. In these expressions, ■Z^i:).'(R)xD-^(R) ~ ('^ 

/v^)'" and 

Zd:^(c) = (7r/2-\/— 5)"^. We finally compare the integrands of the second and the last lines (14-201) 
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and (14-251) of the previous equation that we rewrite as 

ir>j;(R)xDs(ffi) 



X ■ 



Z2RA Jo(n) 



Join) 



dXdY/\\X)K\Yy (^^'+^^'J (4-26) 

DJ(R)xD^(R) 

Zi^^Z^.(K),^.(K)2-m! A(X)A(F(.,)) J^j ^ ^^''^ ' 

Note the sign difference in the two quadratic forms: if X G D^{E?) has eigenvalues Xj, and 
X"" e D^(R) is of the form or then tr {X""^) = -2J2T=i^i = -tr {Xf, and hkewise 

for Y, so that the Gaussian measures match. This justifies a posteriori our change of Z ^ iZ. 

In order to identify the two integrands (the second and fourth hues of (I4-26P ). we notice that 
by definition these integrands belong to L^(R^™) with respect to the measure given by the first 
and third lines of eq. (14-261) . Now we proceed as in [12]: by multiplying F{X,Y) by arbitrary 
polynomials of X, resp. Y, we may multiply the integrands on both sides by arbitrary symmetric 
even polynomials of the Xj or of the Yj. By projecting onto the orthogonal polynomials basis of 
L^(]R^™) with respect to the measure, we deduce that the integrands must be equal. This gives 
the 

Theorem 4.3 For any invariant polynomial function F{., .) and any X"^, Y°- G D"(]R), and X, Y 
the corresponding matrices in D:^{M.), one has: 

I dOe-^'' {x'^OY'^oT) p^^a oY-0'') = (4-27) 

/ 0(n=2m) V ' ; \ / 

0{n=2m+l) 



rG6,n te 



ti 



where 



1 U-' ry ry 

Jac„ L'lRA'^DiiC) 2 2 



(4-28) 



(Jac^)2ZijAZD„^(M)xD„^(M)2"^m! 4™m!Jac^' 

In (I4-27P the dependence on r and the signs ti has been made more explicit, and all ti changed 
into their opposite. 



20 



4.3.3 Examples 

Take as an example the case F{A, B) = 1. Then 



Merging the constant c„ and the triangular integral which decouples into a constant Kn we 
get just a summation over permutations and signs. 



/ 



m 

5Z ^rn[e'^'''''^*''±e-^^''"''^''^] 



ii^n f det2cosh(27X,Fj) 

A(X)A(F) ^ \ det2sinh(27XiFj) 



(4-29) 



with 



Then with 



we obtain 



i^n = c. / e-^*=^(^^^). (4-30) 

^ n 

{/ N 'm{m-l) 
if n = 2m 
^7,™^ (4-31) 
f^j ifn = 2m + l 

n,"ri^(2i)! .r _ ^ 
^n = <l P7ri5' (4-32) 

^ !)■ if n = 2m + 1 . 

This is exactly what was obtained by the Duistermaat-Heckman theorem in sect. 1.3 and serves 
as a check of our formulae. 

5 Correlation functions over the symplectic group 

In this section we repeat the analysis made in section HI in the case related to the symplectic 
group Sp (2m) of 2m x 2m symplectic matrices and to theorem 13. 2[ Following the same steps 
we perform the separation between "angular" and "radial" variables of matrices in both sides of 
eq. fl3^ . 



21 



5.1 Diagonalization of real quaternion antiselfdual matrices 

We consider the set -D^(EI) of real quaternion diagonal m x m matrices whose diagonal elements 
are real quaternions proportional to 62 (see f l2-12p ). with the Lebesgue measure: 

m 

dX:=l[dXi (5-1) 

i 

Appealing again to Cartan's theorem, as we did in sect. 1.4, any real antiselfdual quaternion 
matrix Q G QAm may be written under the form 

Q = SXS^ (5-2) 

where 5 G Sp (2m) and X G D^^{M). 

The Lebesgue measure in QAm is then: 

dQ = Jacf^, A\X) dS dX (5-3) 

with (see appendix ICj) 

"^''''2- = m!n7=i(2j-l)! ^^'^^ 

and A(X) as in (ICTD . 

This decomposition is unique up to a permutation of the m eigenvalues, up to a change of sign 
of each eigenvalue independently, and up to multiplication of 5* by a diagonal quaternion matrix 
V" G Sp (2m) whose diagonal elements Vi satisfy 

Vi = cos 9i + sin 6*462 ; Oi G [0, 2n) . (5-5) 

The latter matrices generate a group isomorphic to O (2)"*. In other words, Q = SXS'' provides 
a 1-to-l mapping between QA^ and Sp (2m) x D^'(m)/{0 (2)"* x 6^ x Z™). 

5.2 Schur decomposition of complex J-antisymmetric matrices 

Let U'^(2m) be the subgroup of U (2m) unitary group satisfying the condition 

U-^ = U^ = JU^r^ . (5-6) 
with the induced normalized Haar measure. We will call these matrices twisted symplectic matrices. 
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Define also Tg*^ to be the set of 2m x 2m strictly upper triangular J-antisymmetric complex 
matrices, with the Lebesgue measure: 

dT := Yl dReTijdlmTij (5-7) 

i<j 
i+j<2m+l 

The 2m x 2m J-antisymmetric complex diagonal matrices of D2m{'^)i (see sect. I4.2p . are also 
J-antisymmetric and come with the Lebesgue measure eq. f l4-8p . Then we prove 

Proposition 5.1 Any J-antisymmetric complex matrix M may always he written as 

M = U{Z + T)U^ (5-8) 

where U G U'^'(2m), T e T/^ and Z G Di^{C). 
The Lebesgue measure in J An is then: 

dM = Jac^^ \A{Z)\'^dU dT dZ (5-9) 
where again A{Z) = Y[i<j{^i " ^j) YIT ^* '^"'^ 

Jacu; = 2-2-Jacf^ (5-10) 
Proof: Similarly to the J-antisymmetric case one can see that 

Det (A - M) = Det (^-A - JAf j) = Det (A + M) (5-11) 

so that the eigenvalues come in pairs (A,— A). One may reorder them to make Z as well as M 
J-antisymmetric and then the constraints on U and T follow. Again the computation of the 
measure follows the lines of (TB]. □ 

This decomposition is unique up to a permutation of the m different eigenvalues, up to changes 
of sign of the m eigenvalues and up to multiplication of f/ by a diagonal matrix V G U'^(2m) whose 
elements are on the unit circle. 

In other words, M = U{Z + T)W provides a 1-to-l mapping between JA^m and U'^(2m) x 

5.3 Symplectic and triangular matrix integrals 
5.3.1 Radial and angular integrals 

Consider the diagonal decomposition of the real antiselfdual quaternion matrices and the Schur 
decomposition of the J-antisymmetric complex matrices. Using them we will rewrite both sides 
of eq. (13-161) . With the measure d/i((5i, Q2) defined in eq. (l3-7p . we have the 
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Theorem 5.1 ^4 matrix integral over QAm x Q-^m can be decomposed into a "radial" and an 
"angular" part using the diagonal decomposition Q = SXS'^ . We have: 

[ dfxiQi,Q2)FiQi,Q2) = (Jacf^)' / dXdYA^iX)A\Y) x 

Jsp{2m) 

Proof: The theorem follows from what is explained in section 15.11 □ 

Notice that one of the two symplectic matrices decouples and so this part of the integral gives 
1. The remaining symplectic integral runs over the relative angular variables. 

Theorem 5.2 A matrix integral over JA2m can be decomposed into a "radial", an "angular" and 
a "triangular" part using the Schur decomposition M = U{Z -\- T)W . We have: 

I dfx{M)F{M,M^) = Jac^^™ /" dZ\A{Z)\^ x 

g-tr(^z2+^z-)g-^tr(z'z) f dT F {Z + T , Z* + T^)e-^ i^''^) . 

Proof: The theorem follows from what is explained in section 15.21 □ 

Notice that the twisted symplectic matrix decouples and so this part of the integral gives 1. 
Only the triangular and radial parts remain. Notice also that the measure for the triangular and 
the radial part factors out and so only a Gaussian measure remains for the triangular part. 

5.3.2 Relating integrals over symplectic and triangular matrices 

Just as in the orthogonal case, we first observe that the integral over the symplectic group 

jSp{2,n).= I d^e-^*'-(^)^^^'F(X,^r^t) (5.12) 

>/Sp(2m) 

with X, y G D^(]HI), is a completely symmetric and even function of the variables Xi and of 
the Yi, i = 1, ■ ■ ■ ,m, since permutation and sign changing matrices may be absorbed into the 
symplectic matrix S. Then the same considerations as in sect. 4.3.2 apply when we want to use 
eq. ( I4-I8p . The Gaussian measure still gets the wrong sign, and the same change of variables 
Z ^ iZ must be used. Then 
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Theorem 5.3 For any polynomial invariant function F(., .), for any X", Y"" G -D^(]HI) and X, Y 
the associated matrices in Do^iM.), one has: 



2m\ 

'A(X)A(F) 



n lib 

X [ dTFitX + T,iY^^,t)+T^)e-^''('^''' 
JtL 



with 



;-l)-JacU-Z2Q^Zn. 



Proof: The proof starts from eq. (13-161) . makes use of theorems 15.11 and 15.21 and then follows the 



same steps as the proof of theorem 14. 3 [ including a change of variables Z iZ, a symmetrization 
in the variables Xi and Yj and the use of orthogonal polynomials. □ 

5.3.3 Examples 

Let's take as an example the case F{A, B) = 1. Then 

_ K)A(Y) 



^Sp(2m) ±^ ™ t_i ^l-^j^l^j 



m T-rm 



K \- \- . TJ , H^kf:^!!^^ K det2sinh(27X,y,: 



with 



Then with 



we get 



K2m = C2m[ e"^*^ (^^^) . (5-14) 

2m 

— 1 (5-15) 

i^'^™ = 2-(™^+2-) ^^"^^^^^ ~ (5-16) 

which reproduces again the Duistermaat-Heckman result. 

In the examples considered in this section and in sect. 3.3.3, the triangular integrals are just 
constants, which is not the case in general. We are going to present an explicit formula to compute 
them. 



25 



6 J/ J-antisymmetric triangular integrals 

In order to compute the correlation functions in the orthogonal and the symplectic group we need 
to compute explicitely various kinds of triangular integrals. For the orthogonal, resp. symplectic, 
case we have to compute integrals over J-, resp. J-, antisymmetric strictly upper triangular 
complex matrices. We shall unify both kinds of integrals into one formalism and explicitely 
perform the integration. 

6.1 Preliminaries to the integration 

The type of integrals we are interested in are of the form 



/ F{X + T,Y + T^)e-^''^^^> dT 
J J 



{6-1] 

I J 

where n is the matrix size, J stands for J or for J, and jj refers to whether we integrate over J 
or J-antisymmetric triangular matrices, and X and Y are J- (or J-) antisymmetric diagonal real 
matrices. Since the measure is Gaussian it is more convenient to normalize the integrals 

(F(X + T,F + Tt))^ = ^-^ ^ — ^ — i- . (6-2) 

■J 

From now on we set 27 = 1, in order to make the propagators simpler. 

The typical functions F[A, B) we want to use are constructed from resolvents and twisted 
resolvents J (;j~4) ■ These functions are not allowed in general by the analytical continuation 
theorems, but this is not a problem if we consider that their x series expansions are generating 
functions of invariant polynomials. An example of such a function is 

T 



tr TT^^J 



1 



tr 7^^^ J ^--^ 7T^^ J 




1 V f 1 Y ^ 1 



X2-{X + Ty \y2-{Y + T^)J \x3-{X + T)J ^y^-{Y + T^))- 

The procedure we use to compute this integral consists in integrating over the last column 
(and by symmetry, over the first row) of the triangular matrices, so as to find a recursion on the 
size n of the matrices, which takes n to n — 2. 

Define the submatrices X and Y by 

X = diag(a,X,-a) ; F = diag(/?, F, -/3), (6-3) 

26 



and the ^7- antisymmetric upper-triangular matrices T of size n — 2 



( 





T 



T 

... / 



V 



T 

n— 1 n 

/ 



(6-4) 



With these definitions and the relations 

oo 



x-(X + T) x-X'^\' x-X 



- Y ^ _ y) ^ ( -r _ 



n=0 



-T 



1 



x-X / x-X 



(6-5) 



we can expand the resolvent of size n in terms of the resolvent of size n — 2 and of the variables 
to be integrated out 



1 



x-{X + T) 



<^i, 1^.7,1 1" ^i,n^_i,n [ 1" ^i,l^j,n ^1 



X — a X + a 



+ (1 - <^i,l - <^i,n)(l - <^i,l - Sj^n) 



X — a X + a 
1 



a; - (X + T) 



-|-(5i,l(l - - 5j^n] 



1 



-f (1 - - 5i^n)Sj,n 



X — a 



1 



.=2 VX-(X + T) 

1 



1 



X + a 



1 



X — a 



1 



X + a 



. (6-6) 



2<fc</<n-l \^ ~ ~^ '^^ J k,l 

Notice that Ti „ in the J-antisymmetric case is identically zero, a fact that will be accounted for in 
the following. In both cases, the only independent integration variables are the matrix elements of 
the first row. Their propagators are read off the Gaussian weight, which is, in the J-antisymmetric 
case 

Q-hU^i{\Ti,.\^+m,n\^) = Q-u^i i^L.P (6-7) 

while for the J-antisymmetric case it is 



The independent nonzero propagators are thus 



(6-8) 



J J , 2<iJ<n-l 
1 + b 



(6-9) 
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where b = ~1 (resp. 6 = +1) for the J-antisjTumetric (resp. J-antisymmetric) case, so that in 
the b = —1 case, the propagator for Ti „ is zero, as it should. The other propagators encountered 
in the integration result from the symmetry properties 

'T^,nTl^) = {Ti,TI) = -J., {T^,nTl) = - . (6-10) 

This is what is needed to perform the first step in the recursive computation of the triangular 
integrals. Let's take the simplest mixed case. 

6.2 Example: Morozov-like formula 

The simplest case involves two resolvents. Define the two functions 

Fi''\x,y,A,B) = ti ( rrl+l 

+ ^ ' ^ \x-{X + T) y-{Y + T^)J 



T 



The second one, is twisted by the action of J . Using that, for A and B two ^/-antisymmetric 
matrices (such as X, Y , T or T), 



X — A J X + A 

1 V. =. 1 



i.e. J7 (^3^)"^ JT" = , one sees that F^\x,y, A, B) = bF'i\x, —y, A, B), thus it suffices to 
carry out the integration over the last column and first row of only. 
The computation goes as follows 

\ r 11 111 

l)(n-2) 



' x-(X + r) y-(F + Tt) ' +^ 



in) 



11 11 

1 + t; + 



x—ay—P x+ay+P 

1 1 



^x^ -a^y'- + \ _ (X + f ) y - (r + ft) 

I I I ( 1 1 



(n-2) 



X - ay - P \\x - {X + T) y - {¥ + T^) 



(n-2) 



+ ..^. ( I ^^-^ ^-^1 > (T.nTl) (6-12) 



x + ay + p \ \y-(Y + T^)x-{X + T) 



1 I 1 \ / 1 



(n-2) 



■> ^ [n 2) 
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Inserting the propagators given above we find 
1 1 



tr 



x-iX + T)y-{Y + T^) 



1 



1 



+ 



X — ay — P 
1 1 



1 



(n) 
1 



a; + a y + /3 
1 



tr 



tr 



-J 



1 



a; - (X + T) 1/ - (r + Tt) 

T 



1 



1 



x-{X^T) \y-{Y + T^) 



J +b 



(n-2) 

.(6-13) 



(n-2) 



We have spht the two terms in the ^Ti „tJ^ = 1 + 6 propagator in the following way: the weight 
1 goes together with the untwisted minimal cycle and the weight b with a twisted one. This is a 
general rule as we shall see later. 

Notice that we need both functions F± in order to close the recursion relation. Defining the 
column vector V„ = -FI"'*)^, we obtain a recursion formula for the two functions in the form 



Vn = M{x,y,a,p)Vn-2 



(6-14) 



where 

M{x,y,a,(3) 



1 + — — I (1 + ^^ 

x—a y—f) I \ x+a y+P 

1 1 



This structure will appear in the general case. 



1 1 

x'^—oi^ y'^—fj'^ 
1 I \ (l ^ I 1 



x~a {-y)-f3 



x+a {-y)+P 



(6-15) 



6.3 Last-row/first-column integration: general case 

In the general case, the recursion involve combinations of correlation functions conveniently la- 
belled by graphs. 

6.3.1 Basis of correlation functions 

In the example above, we had to mix correlation functions with two resolvents with correlation 
functions with a lesser number of resolvents, in order to write recursion relations. This will still 
be necessary in the general case, and we shall see that the basis of correlation functions that we 
have to consider is conveniently labelled by tetrads uj = {a, r, s,t} E &r x Gf{ x x Z^, made 
of two permutations of Qn of R objects and two sets of R signs. This integer R will turn out to 
be the maximal number of resolvents of X-type appearing in the correlation function. 
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To each such tetrad, we first associate an oriented bicolored graph G in the following way: G 
has 2R vertices, R of each color; the zth black (resp. white) vertex carry a sign s{i) (resp. t{i)), 
i = 1, ■ ■ ■ , R. An oriented edge connects each ith black (resp. white) vertex to the a{i)th white 
(resp. r~^(j)th black) vertex. The graph is thus made of bicolored cycles. 
According to these rules, the graph representation of the tetrad 

u = {(13)(24), (1)(243), (+, +, +, -), (+, -, -, +)} (6-16) 




For our purposes we are only interested in the relative signs between vertices in each cycle. We 
say that two tetrads u and u' are equivalent if they are equal up to independent global signs in 
each cycle. We take as representatives of the equivalence classes the tetrads with s{i) = + for the 
black vertex with the smallest i in each cycle. 

There exists also a graphical representation for the equivalence classes. We call diagrams these 
new objects. Since we only care about the relative signs between vertices in the same cycle, we 
replace the Z2 variables by symbols representing changes of sign 0. We represent this with a short 
bar across the edges indicating the change of orientation of the edge. In the diagrams, the vertices 
will be called dots and the edges will be called links. The equivalence class 



M = [{(13)(24),(1)(243), (+,+,+ 



) / ) V I ) ) ) 



+)}] 



(6-17) 



is thus represented by the diagram 



4 
rOi 



^This operation is accompanied by a change of orientation of edges, and the cycles are no longer oriented. But 
the original class of oriented graphs may be reconstructed from these diagrams. 
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In appendix [E] we prove that the space of equivalence classes of tetrads is isomorphic to the space of 
permutations of size 2i?, so we can label the equivalence classes [uj] by an associated permutation 
^ ^ ©2_R- This very important bijection will allow us later to relate the O (n) and the Sp (2m) 
integrals with the f/„ group integrals studied in |12j . 

The correlation functions we need turn out to be labelled by these equivalence classes, or 
equivalently by permutations of size 2R. We define in general the [u] component of the basis of 
correlation functions as follows. Consider [uj] = [{a, r, s, t}] and consider two sets of variables, 
{xi, . . . ,xr} for X-type resolvents and {yi, . . . for F-type resolvents. Then we define the 
following function, 

n («(.")^«.. ^ - (n (^) ) j 

for the J-antisymmetric integrals and 



for the J-antisymmetric ones, where p is the number of cycles of the permutation ar"^, and Rk 
is the length of the ki\i cycle. 

The permutations a and r G &rn yield the ordering of the labels 

(y{ik,i) = jk,i , and T~^{jk,i) = ik,i+i , with = ik,r < ik,r r = 2, Rk, (6-20) 

and the signs s, satisfying the constraints 

s(zfc,i) = +1, (6-21) 

together with the signs t define the functions 

x^y _ i if s{ik,i) = t{jk,l) 



Vk,i 



and the operations 



C{jk,i) 



1 if s{ik,i) = -t{jk,i) 

iit{jk,i) = s{ik,i+i) 

1 if t{jk,i) = -s{jk,i+i] 



T(ranspose) if s{ik^i) = —I 

J(dentity) if s{ik,i) = 1 

T(ranspose) if t{jk,i) = —1 

/(dentity) if t{jk,i) = 1 
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that perform the twisting of resolvents. 

The structure of these functions is easily understood from the diagrams associated to vr (equiv- 
alently, [uj])- Each cycle in the diagram represents by a trace; to each dot is attached a resolvent 
if the dot is traversed clockwise by an arrow, and a transposed resolvent if it is counterclockwise; 
finally each change of orientation in the links corresponds to a J or a J matrix. 

The functions are invariant under an independent global twist inside each trace, so the claim 
that our prescription depends only on the equivalence classes defined above is justified. 

Finally the terms 5_Rj.,i and t(j)5_R^,,i in the definition of the functions are the analogues of the 
1 and b appearing in F± in the example of sect. 16. 2[ Here too they come only with the traces 
containing two resolvents, i.e. = 1. 

As an example the [{(1)(2)(3), (1)(23), (+, +, — ), (— , — , +)}] component of the basis for the 
orthogonal case would be 

T 




which is represented by the diagram 



-1 + 




There is a unified representation for the two bases of correlation functions corresponding to 
the orthogonal and symplectic cases. Define Nj = YlkiivtT^ + vtT^) (which is the total number 
of JT" matrices appearing in the traces), and n(s) = Ylf ■s(^) for s G Z^. Then we havt^ 



ltr(^-)n(,n(,n(...,.tr 

itr (:r^-)n(s)n(t)(F5,,({x}2«,M2^,A,i?))^ . 



^t{jk,i)yj,, - B 



^Here we have used eq. (|6-lip in order to remove the J and introduce the signs. 



(6-22) 
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where n G &2r is the associated permutation following appendix [El e, {x}2r and {y}2R are 
also defined in that appendix, and F^^, is the basis of correlation functions found in [Hj for the 
unitary case. The sign can be computed through the limit A, B ^ oo. We prove the last equality 
in appendix [El 

6.3.2 Recursion relation 

Using this basis we have the following theorem 

Theorem 6.1 The functions defined in eq. ^6-18\) and eq. ^6-19\) or equivalently in eq. ^6-2^) 
satisfy the recursion relation 

{<T',T',s\t'} 

where a and (3 are the first eigenvalues of X and Y respectively, and A is the suhmatrix of size 
n — 2 resulting from erasing the first and the last rows and columns of A, and where 



R 



'i)xi + at{a{i))ya(i) + [3' 



-I— r 1 1 

n(^^«,.'W^^«.'W%»),t'(^W) + s{^x.-at{r{^)y^^^-p\ 

Proof: The proof of this theorem is given in appendix [Dl □ 

The first thing to notice is that (see appendix [El for a proof) the matrix M. is again closely 
related to the corresponding recursion matrix M. found in [12] for the unitary case. Recall that 

M^"^X{x}{y},i,r^) = n + /_ . , , (6-24) 

then taking into account the bijection defined in appendix [El we find that 

M^""^ {y}R, «, P) = -^S^ i{^}2R, {y}2R, -a, -(3) (6-25) 

which again connects the orthogonal and symplectic cases with the unitary case. The precise 
relation between the arguments of the two sides of this equation is defined in appendix [El Trivial 
consequences of this fact are the commutativity property of M. 

7W^''\{x},M,a,/3),7U^''\{x},M,e,r/)] =0, (6-26) 
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and the symmetry A4 = Ai^ . 

The recursion relation we just found is vahd for any value of n such that n > 3. The special 
cases n = 2 and n = 1 correspond to the initial condition for the recursion relation in the even n 
and odd n case respectively. 

6.3.3 Initial conditions 

Let us consider first the case n even. For any n > 2 even, the recursion relation is valid. The last 
step for n = 2 requires a slightly more careful analysis. In a 2 x 2 strictly upper triangular matrix, 
the only term is T12. Since the "last column/first row" integration reduces to that of T12, the 
procedure explained in appendix [D] is still valid, and the recursion relation can be naively applied 
just by considering that when n = 0, all the traces are equal to zero in the correlation functions 
(or equivalently taking the strict limit where Xi, yi ^ 00 ioi i = 1, R). This gives us the initial 
condition vector 

{I^){a,r,s,t} = tr {J""-^) U{s)U{t)S„,, . (6-27) 
Explicitely, the two cases J = J and J = J are 

{U){.,r,s,t} = n(s)n(t)5^,. (6-28) 

The case n = 1 corresponds to the case where there is no triangular matrix at all, thus no 

integration. The initial condition corresponds in this case to the vector 

P f R, \ 
{I(){.,r,s,}{{x}. {y}) = tr ( J^^) n t{j,,)6n,,, + n • (6-29) 

k=l\ i^^^ik.iVjk.i J 

7 Correlation functions over O (n) and Sp(2m): the final 
expression 

In this section we find a determinantal formula for the correlation functions in sections H] and [51 
We use need the matrix determinant Mdet defined as, 

n 

Mdet(M)= 5^ (-1)'^J]M,,.(,) (7-1) 

o-G6(n) i=l 

where each Mij is a matrix. This means that Mdet(M) is itself a matrix. The ordering of matrices 
in the product does not matter if the matrices commute with one another. 
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• The O (2m) case 

Take theorem 14.31 for n = 2m and the recursion relation found in section [Hi After some algebra 
we find 

Mdet (e2T^^^^A7^''\{a;}, {y}, tX,, tYj) + e-2^^^^^A7^''\{a;}, {y}, zX^, -zY,)) 

_ ^ / k,j=l,...,m jj 



det (2 cosh i2^XkYj))^^^^^ 
Mdet [e^i'^^''^^M^^\{x},{y},iXk,i t^Y^ 



j=i,...,m 



= V ^ k,j=l,...,m jj s 

det(2cosh(27X,F,)),^^.^,_ „^ ° ^ ' 

where F^^^ is a vector with (2i?) ! components, Mdet is a (2i?) ! x (2i?) ! matrix and J^^ is the vector 
defined in eq. ( l6-28p . 

• The O (2m + 1) case 
Take now the n = 2m + 1 part of theorem I4.3[ With the results of section [6] we find 

/o(2^+i) dOFi^\X, OYO^)e--- (^^^^") 

Mdet (e'^''-^^M^''\{x}, {y}, zX,, lY^) - e-''^^^''^M^''\{x} , {y}, iX,, -lY^)) 

V / fc,7 = l 



j = l,...,m 



det (2sinh(27Xfer,)), .^^^ 
Mdet (e2'>*^^^^^AT^^^({x}, {y}, iX^, i t^Yj] 



tijVii / det(2smh(27XfcF,)),_^.^^_ 

where // is given in eq. (16-290 . 
• The Sp (2m) case 

Finally take theorem 15.31 and the corresponding part of section [61 After some algebra 
Is,,,m,dSF^^HX,SYS^)e~--i^'^'') 

Mdet fe2^^^^^AY^''\{x}, {y},tXk,tY,) - e'^^^^-^^ AY^^'^Ix}, {y},tXk, -lY, 



k,j=l,...,m jJ 



det(2sinh(27Xfcr,)),^^.^,_ 
/ m \ Mdet (e2^*^^'=^:'JW^''\{x}, {y}, zX,, zt.F,)) 
^5 IG'V det(2smh(27X.r,)),,.,..,. '^""""^ ^'"'^ 

where /q is given in eq. (16-281) . 
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8 Concluding remarks 



A remark on contour deformation in matrix integrals 

We first return to the transformation of our original integrals over two real antisymmetric, resp. 
two antiselfdual real quaternionic, matrices into integrals over J-, resp. J-, antisymmetric complex 
matrices. Define the following measure on M„(C) x M„(C), the set of pairs of two complex n x n 
matrices, 

^-tr (i^A/^l^M|+,M,M.)dMidM2 . (8-1) 

We will consider two hyperplanes of M„(C) x M„(C) and the measure on these hyperplanes induced 
by the measure above. 

The first hyperplane is defined by the equations 



Mi = M* = Ai 
Ai = Ai 



1,2 (8-2) 



endowed with the induced measure eq. (13- ip reproduces the two real antisymmetric matrix inte- 
gral of 13.1. 1[ The second hyperplane, which describes one complex J-antisymmetric matrix, i.e. 
JAn{C), is defined by the equations 

Ml = Ml = M , . 

JM = -M^J 

with the induced measure eq. (13-40 . 

In this construction, the real antisymmetric two matrix integrals and the complex J-antisymmetric 
matrix integrals are nothing but the same integral on different hyperplanes of M„(C) x M„(C). 
By the counting done above, these two hyperplanes have the same dimension, and a plausible 
interpretation of theorem 13.11 is that it results from a contour deformation taking the first set of 
matrix integrals into the second one. 

Similarly we can consider two different hyperplanes in the space of pairs of quaternionic ma- 
trices, namely the one defining QAm x QAm 



Mi = M* = Q, 

Qi Qi 



1,2 (8-4) 



with the induced measure eq. (13-71) . and the one defining JA2m 

(8-5) 



Ml = Ml = M 



JM = -M^J 
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with the induced measure eq. (13-131) . Both hyperplanes have the same dimension and we may 
again interpret Theorem 13.21 as resulting from a contour deformation. 

This apparently much simpler and intuitive approach has the drawback of neglecting conver- 
gence issues. This is the reason why we chose to prove our results by use of loop equations. 

Other choices of deformations to other hyperplanes are also conceivable. 

Comparison with the Duistermaat-Heckman form 

First observe that we may rewrite the main theorems 14.31 and 15.31 together with theorem 4-1 in 
[T2] into a unified form involving an integration over a set Tq of complex triangular matrices, 
as already mentionned in the Introduction. When G is 0(n), Sp{2m) or U{n), Tq corresponds 
respectively to the set of J, J-antisymmetric, or unconstrained, strictly upper triangular matrices. 
Moreover, we notice that in each case, the set Tq is precisely the derived ideal [b, b] =: n+ of the 
Borel subalgebra b associated with the choice of Cartan algebra made above. This is in fact the 
space generated by the positive roots n+ = ©a>o 0°. 

It is thus natural to conjecture that an analogous formula holds for any compact group, with 
the identification of Tq with the subalgebra n+: 

Conjecture 8.1 For any compact group G 

[ dfiF(X^fiF<^^]-l)e-*^(^"^^''^") 
Jg 

+ tr{Xw{Y)) r 

This form of our result has to be confronted with the form given by Duistermaat-Heckman's 
localization theorem, eq. (12-41) . Note that the integration over n"'' plays here the role played in 
sect. 12.21 by the integration over the "fluctuations" A G in ( 12-4^ . This points to a possible 
much more compact and geometric derivation of our results. 

Other comments 

First, it is remarkable that the recursion on n for triangular integrals involves the same matrix 
Ai for all cases U(n), 0(n), Sp(2n). Only the initial conditions differ. This fact needs to be 
understood, and it shows that the matrix A4 is universal. Moreover, its commutation properties 
suggest the existence of some underlying integrable structure. The symmetries of the group 
under consideration are reflected in the symmetries of the spectral parameters at which A4 is 
evaluated. The initial conditions also seem to have such symmetries, and it is remarkable that 
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those symmetries are reminiscent of root lattices of size R or 2R (we started with a root lattice of 
size n or 2n) . This suggests a duality between R and n, similar to that of the "supersymmetric" 
method of evaluation of determinantal correlation functions [20]. Remarkably, the triangular 
matrix ensembles as those considered here seem to play an important role in generalizations 
of the so-called Razumov-Stroganov conjecture. Indeed multidegrees of the corresponding matrix 
varieties are solutions of the quantum Knizhnik-Zamolodchikov equation based on the root systems 
of type A, B, C and D [21], thus pointing again towards some possible integrable structure. It 
might be also interesting to attack the "angular" integrals considered in this paper with character 
expansion techniques, see [22] for recent references. 

Our last comment is that it would be highly desirable to know how to compute integrals like 
(11-11) on other orbits. For example, little is known about the integral over the 0(n) group when 
X and Y are symmetric real matrices (see however [23]). 
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Note Added: 

At the date of this resubmission, conjecture 18-61 has been recently proved in [23] giving a full 
meaning to each of the elements appearing in the conjecture. The role of the borel subalgebra 
and the complexification procedure is thus quite well understood. 

Appendices 



38 



A Quaternions 



We review here some well known facts about quaternions to fix our notations, which follow basically 
those of [16]. We shall only consider the set of real quaternions, which is the algebra over M 
generated by 4 elements: the neutral element cq, which by an abuse of notation we often write 
Co = 1, and Cj, z = 1, 2, 3, 

q = g(o)eo + g(^)ei + g^'^ea + q^^^e^ q^"^ E M (A-1) 

with multiplication ef = 616263 = — 1, from which it follows that 6162 = —6261 = 63 and its cyclic 
permutations. One may represent 60 by Id2 the 2x2 identity matrix, and the 6j in terms of 2 x 2 
Pauli matrices by 

Ci = -iai . (A-2) 

The conjugate quaternion of q is defined as 

g = g(0)l-g(l)ei-g(2)e2-g(=^)63 . 

(This is also called hermitian conjugate, which is justified by the fact that Pauli matrices are 
hermitian). Note that gg:= = the square norm of the quaternion, 

and hence q ^ ^ has an inverse q~^ = q/lql"^- Real quaternions form a non-commutative field. 
Note also that conjugation reverses the order of factors of a product (^152) = 92^1 • 

Quaternionic matrices 

We now consider matrices Q with real quaternionic elements Qij, i, j = 1, ■ ■ ■ ,m. Alternatively, 
using (lA-2p . one may regard also Q as a 2mx2m matrix with 2x2 blocks made of real combinations 
of Id2 and the Pauli matrices. One may apply to Q the same conjugation as defined above. One 
may also transpose Q. The dual of a quaternionic matrix Q is the matrix 

iQ%j = QJ^ . (A-3) 

(This is also the hermitian conjugate of Q in the usual sense.) A real quaternionic matrix is 
thus self-dual if 

Q^ = Q = Q^ = (Q,j) = (g,,) . (A-4) 
A real quaternionic matrix is anti self-dual if 

Q^ = Q^ = -Q; (A-5) 
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it is thus anti-hermitian. In particular, its diagonal matrix elements are such that Ql^ = 0. 
On quaternionic matrices, we may define the ordinary trace 

m 

tr {Q)=J2Q^^ ' 
i=l 

which is in general a quaternion, or 

m 

tro (Q) = 2 ^ Qi°) = tr (Q) + MQ) (A-6) 

which is a scalar. Note that tro [Q) is nothing else than the trace of the corresponding 2m x 2m 
matrix. 

Symplectic group Sp (2m) 
Let H be the space of real quaternions. Consider the hermitian form on H*" 

m 

{x, y) = ^ XiVi . (A-7) 

i=l 

The compact unitary symplectic group Sp(2m) is defined as the invariance group of that form, 
and is thus the group of m x m real quaternionic matrices Q such that 

Q^Q = I or g^g = /. (A-8) 

These matrices may be called unitary real quaternionic matrices. The Lie algebra of Sp(2m) is 
generated by real quaternionic matrices A satisfying the infinitesimal version of (A-8), 

+ A = Q or A^ = A^ = -A , (A-9) 

hence by antiselfdual real quaternionic matrices. 

B Loop equations I 

In this appendix we show how to compute loop equations in the real antisymmetric two matrix 
integral and in the J-antisymmetric complex matrix integral. 
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B.l Loop equations for the 2 real antisymmetric matrix integral 
B.1.1 Loop equations 

Schwinger-Dyson equations, also called loop equations in the case of matrix integrals, merely 
amount to saying that the integral of a total derivative vanishes: 

= J2[ dAdA, A,).,e-(^^?+^^i+-^^^^)) (B-1) 

i<j 

where f{Ai,A2) = —f^{Ai, A2) is any sufficiently regular matrix valued function; in particular / 
can be any non-commutative polynomial in Ai and A2, and may contain also product of traces of 
polynomials. 

The loop equation thus turns into an equality between expectation values: 

= (tr ((aiAi+7A2)/(Ai,A2))) (B-2) 

where 

Notice that Ki{f) is linear and satisfies Leibniz rule: 

K,ifg) = K,{f)g + fK^{g) . (B-4) 

The most general / we shall consider is of the type: 

1 ^ 
f{A^, A2) = -mA^, A^) - Fo{A,, A^f) J] WA,, A^)) (B-5) 

r=l 

where Fq is an odd degree non commutative monomial of Ai and A2: 

Fo(Ai, A2) = 4°'°^?°''^2°''^?°'' • • • a'°'''"-'a';"'''A'°''° (B-6) 
and each Fr with r > 1 is an even degree non commutative monomial: 

and we call deg(/) the total number of matrices Ai + the total number of matrices A2. 
Then compute: 

po ko,q-l 

K,{Fo -F^) = E E E ■ ■ ■ A'r ArUA>"^~' ■ ■ ■ A'^n (B-8) 

5=1 m=0 i<j 

( aIo.o 4^0,1 4^0,1 ( Ako.q-m-l .«o,q A^O:Po\ 
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( aIo,0 Ako,i aIo.i a'-i-'^ Am\ ( A^'0,q"rn~l .lo,q 

~\^2 ^1 ^2 •••^2 ^1 ^2 • • • 



I / 4^0,0 /lfco,l /iio,l /iniN / /|fco,<7-m-l ««o,q /I' 

T"l^2 ^1 ^2 • • • ^2 ^1 ^2 • • • ^ 



2 Jjj 



EE 

(j=l m=0 



J] J] |tr [A^AT^^ArA, 



lo,a 4^0,1 4«o,i /,fco,2 4^0,9-1 ( Ako,q-m~l j^o^q ^ 



^2"''-^™ tr [Al 



-(-1) 



ko,i-\ hfeo.q-i+m+'o.oH \-lo,q-i 



9-1^^ / ^rrt^'o,q-i ^^0,2 ^10,1 ^to.i ^'0,0 



fco,2 /l'o,l /l'=0,l /I'O.O Ako,q-m-l ,lo 



j^'-^,q 



This equality is known as the spht rule. 

Then we have for any antisymmetric matrix C: 

Pr kr.q — 1 



( A^"^'^ A^^-^ A^^''^ A^r.q~i Ain\ 1 Akr.q—m.~l 

-y/i-l /I2 /li • • • ^2 ^1 )sj ^ij 1^1 

\ ^ \ ^ ( J^r,l J^r,l J^r,q-1 ^ J^r,q—'m-l j^r,q J^' 

/ _j / _j \121 211 2 2 



5=1 m=0 i<j s 



Pr kr.q — 1 



A ''r,Pr \ 

^2 Jjs 



A'-r.q A'-r.Pr\ 

J-I2 ■ ■ ■ ^2 J is 



(B-9) 



q=l m=0 

This equality is known as the merge rule. 



Due to Leibniz rule and using repeatedly the split and merge rules, we find that if / has the 
form of eq. (lB-5p then Ki[f) is a linear combination of monomial invariant functions of degree 
< deg(/) - 1. 

The loop equations read: 



(tr (aiAi/ + 7A2/)) = i^i(/) 
(tr {aM + ^A^f)) = K2U) 



(B-10) 



or equivalently: 



(tr {A,f)) 



(tr iA,f)) = ^K,{f)-lK,{f) 



(B-11) 



B.1.2 Polynomial invariant functions 



F{Ai,A2) is a monomial invariant function of two antisymmetric matrices Ai,A2, if it is either: 

(B-12) 



F = 1 

or F{A,,A2) = tr (^1/(^1,^2)) 
or F{AuA2) = tr (^2/(^1,^2)) 
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where / is of the following form: 

R 

f{Ai, A2) = Fo{A,, A2) n (FriA,, A2)) (B-13) 

r=l 

where Fq is an odd degree non commutative monomial of Ai and A2: 

Fo(Ai, A2) = • • • a'°'''"-'A'1"'''A'^''° (B-14) 

and each F^ with r > 1 is an even degree non commutative monomial: 

Fr{Ai, A2) = J^{-^ J^^-^ J^{'^ . . . A^--^--^ a\^'^^ A^;-"^ (B-15) 

and we call deg(F) the total number of matrices Ai + the total number of matrices A2. Notice 
also that / can be antisymmetrized without changing F, and thus / can be taken of the form of 
eq.dEIS]). 

Notice that if deg(F) is odd, we have: 

(F) = (B-16) 

If F = 1, i.e. if deg(F) = we have: 

(1) = 1 (B-17) 
and if deg(F) > 0, and F = tr (Ai/), the loop equations imply: 

{F) = ^K\{f) - lK2{f) (B-18) 

where the right hand side is the expectation value of a polynomial invariant function of degree 
< deg(F) — 2. And if deg(F) > 0, and F = tr (^2/), the loop equations imply: 

(F) = ^K2{f) - ]K,{f) (B-19) 

where again the right hand side is the expectation value of a polynomial invariant function of 
degree < deg(F) — 2. 

In other words, the loop equations allow to compute every expectation of polynomial invariant 
functions by recursion on the degree. 

Notice also that the expectation value of any monic monomial invariant function is a polynomial 



in a, ^ and |. 
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B.2 Loop equations for the complex J- antisymmetric matrix integral 
B.2.1 Loop equations 

Similarly to the previous section, loop equations, in the case of a complex J-antisymmetric matrix 
integral, can be written: 

(B-20) 

where /(M, M^) = —Jf{M, M^^)^ J is any sufficiently regular matrix valued function, in particular 
/ can be any non-commutative polynomial in M and M\ and may contain also product of traces 
of polynomials. 

The loop equation thus turns into an equality between expectation values: 

{K,{f)) = (tr ((aiM + 7Mt)/(M,Mt))) (B-21) 

where 

Notice that Ki{f) is linear and satisfies Leibniz rule: 

K^{fg) = K,{f)g + fK,{g) (B-23) 

The most general / we shall consider is of the type: 

1 ^ 
/(M, Aft) = -(Fo(M, Mt) - JFo(M, M^f J) J] tr (F,,(M, Aft)) (B-24) 

where Fq[M, M"^) is an odd degree non commutative monomial of M and M^: 

Fo(Ai, A2) = Mt'°-«M^o-Wt'«'iAf'^»>2 _ ^ ^ Mt'o-fo-iM^^o-foM^'O'TO (B-25) 

and each F^(Af, Af"^) with r > 1 is an even degree non commutative monomial: 

FriM, Aft) = M^^'-' M^^^'' M'''-'^ . . . Mt^'-.f'-iM^'-'^'-Mt'-'''- (B-26) 

and we call deg(/) the total number of matrices M + the total number of matrices Aft. 
We have 

K,{Fo-JF^J) 
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g=l m=0 i<j 



i,j V ■ ■ ■ / n+l—i,n+l—j 

\ ' j,i V J n+l—j,n+l—i 

+ ( M^''°-° M'"^-^ . . . M^''"-'"-^ M'^) (m''°'''~'^~^M^''°-'' . . . M''°'™ M"^'"-™ ) 

V ■ ■ ■ / j,j \ ' ' ' / n+l—i,n+l—i 

PO ko,q-l 

= [tr (M"^'0'«M^"'^ . . M^'^'-'-W™) tr (M^"'«^™"^M^'°'« . . . M^^'TOM'^'O'to) (B-27) 

q=l m=0 

_(_]_y(),oH hio,9-i+A:(),iH hfco,q-i+m^^ ^ j^m j^t«0,q-l J^ko,i j^\lo,o ]^ko,q-m-l j^\lo,q ]\^'=0,po ]\^t'o,po j j 

Notice that the spht rule for J-antisymetric complex matrices is identical to the split rule for real 
antisymmetric matrices eq. flB-Sp . 

Similarly, we compute the merge rule (where C = —JC^J is any J-antisymmetric complex 
matrix) : 

K,{CtT (Fr)) 

Pr ko.q—1 

q=l m=0 i<j s 

Pr fco.q — 1 

= ^ ^ tr (M'^'-'I . . . M^''-'''-W" CM'''-'9"'""^M^"''-'^ . . . M^'-'P- M'^^'-'P'-) (B-28) 

q=l m=0 

And again the merge rule for J-antisymetric complex matrices is identical to the merge rule for 
real antisymmetric matrices eq. flB-9p . 

We conclude that the expectation values of invariant polynomials of M and are entirely 
determined by the same recursion relations (on the degree) as the expectation values of invariant 
polynomials of two real antisymmetric matrices. This completes the proof of Theorem I3.1[ 

B.3 Symplectic case 

The procedure to obtain the loop equations for the real quaternionic antiselfdual two-matrix inte- 
gral and for the J-antisymmetric complex matrix integral and to prove Theorem 13.21 is completely 
analogous to the one above. 
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C Calculation of Jacobians. 



In this appendix we are going to detail the main steps for the computation of the Jacobians fl4-3p . 
f l4-12p . fl5-4p and f l5-9p . For this purpose we will need one of the special limiting cases of the 
Selberg integral called the Laguerre limit (see for example |16j ) 



vl<i<jr'<n / \fe=l 

n-1 



-|-|- r(l + 7 + J7)r(a + j7) 

r(i + 7) 



The two values of this integral we need are 



^(il'-)=-!i^ilW (C-2) 
/(^,l,m) =m!^^f[(2j-l)! (C-3) 



o 



The procedure is essentially the same for the four cases. Let us show in detail the first one, Jac! 
The following relation holds true by the block- diagonal decomposition shown in section 14.11 



/ dAe^K-) = jac^ r n (^?-^?)'{ 



.Tr(^) _ T^^o r TT t^2_^2^2j nr=ie ^■'^dxfc ifn = 2m 

YYk=i xle'^i'dxk if n = 2m + 1 



l<i<j<ni 



- Jac„ j ^^1^^^^^ ifn = 2m + l. ^^'^^ 
Computing the Gaussian integral on the left hand side we find 

jac„ - [V^) \ if n = 2m + 1 

which gives exactly the expression in equation l4-3[ 

The Jacobian Jacg^ is computed with the same technique from a real quaternionic antiselfdual 
Gaussian integral, 

i) iV^r^'-^'^ = I dge"" (-) = Jacf^J (I l,m) . (C-5) 

In order to compute the two remaining Jacobians Jac^"^ and Jac^^ we need to introduce 
two new matrix ensembles. Consider the Hermitean J / J-antisymmetric one-matrix model. By 
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the J/ J-antisymmetry of these matrices we know that they can be put into a triangular form 
by a twisted orthogonal or symplectic matrix respectively. By hermiticity we conclude that the 
triangular part of this Schur form will be zero, and the diagonal part (eigenvalues) is real. It is 
easy to argue that the Jacobians for these transformations have to be the same as the ones we 
seek in sections 14.21 and 15.21 This allows us to write the following 



dHe 



. n(n — 1) _ 

(v^)^™''^™2"'if = J 



/ (|, 1, m) ii J = J and n = 2m 
Jac^"" <j / (|, 1, m) if = J and n = 2m + 1 (C-6) 
2^™/ (|, 1, m) if J' = J and n = 2m 



which gives eq. fl4-12p and eq. fl5-10p . 

D Proof of theorem 16.1 



In this appendix we use the graphical representation of the basis of correlation functions introduced 
in section [B.3.1l to prove theorem l6.1[ The idea is to identify all possible occurrences of elements of 
the first row and last column of the T matrix (and vice versa for T'^) by means of the decomposition 
of eq. (16-61) : then to use the constraints coming from (i) the triangular structure of these matrices, 
(ii) the contractions of indices within traces, (iii) the propagators eq. (16-91) and eq. (16-lOp . to 
represent the result of the integration in a graphical way, leading to the recursion formulae. 
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D.l Last column/first row integration 



Take the functions defined in 16.3.11 and their graphical representation. We first rewrite eq. (16-61) 
in a shghtly reshuffled form 



x-iX + T) 



X — a ' X + a 
1 



X — a 



X + a 



(D-1) 




\a:-(X + T). 

where a and (5 are the first eigenvalues of X and Y respectively. We substitute this expression 
for each resolvent in the integrand, and perform all possible "contractions" of the Ti fc, Ti ,„ and 
T^^n variables by means of the propagators (16-91) . This can be represented as operations on the 
diagram associated to the given function. Note that the terms in the last four lines of (ID- II) still 
contain a resolvent (of size n — 2), while those on the first two lines do not. Let us now enumerate 
the operations corresponding to each term in the expansion eq. (ID-ll) : 

• Operation 1: The term on the first line, which singles out one Ti„ variable, removes one 
resolvent from the integrand, which is represented by erasing a dot in the diagram. Since 
Ti „ can only be contracted with the appearance of this term forces the erasing of a dot 
of the opposite color, by another application of Operation 1 on a y-type resolvent, somewhere 
in the diagram. Since this Ti„ appears in a trace, its left and right neighbouring resolvents 
must have a with one matching index 1 or n. The operation of erasing dots leaves pairs 
of free links with only one dot at their end, carrying such a T"^ variable; their role and their 
weight will be reconsidered in Operation 3. The same applies to the other erased dot. Let 
us now perform the contraction of the selected Ti^nT^i pair, giving a factor (1-1-6). The 
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graphical representation is, 

( ^ ' ^\ 




1 1 



( 



\ 



These pairings have to be performed in all inequivalent ways. 

Operation 2: The two terms in the second line play a similar role. They also remove a 
resolvent, which is again represented by erasing a dot. This forces one of the neighbors to 
be replaced by a similar term. This will be represented by the operation of erasing a link 
and its two adjacent dots. The possible configurations are 




Xi - ayj - P 

1 1 

Xi + ayj + (3 

1 1 

Xi-ayj + P 

1 1 



Xi + ayj - (3 

The signs in the last two equations come from the twist of the erased link. This operation 
also leaves some free links. 

Operation 3: All the remaining terms do not remove resolvents so they do not erase any dot 
in the diagram. Instead, they represent cuts in the links, since each Ti ^ and „ forces a 
T^i or tIi^ in a neighbor. We must also consider here all free links created by erasing dots 
in Operations 1 and 2. As discussed above, these terms contain also a T or a T"!" variable at 
the end of the free link, and will contribute to the weight. Graphically we have 





for the cutting, and 





for the free hnks coming from Operations 1 and 2. In the right column, the bar across the 
free hnk indicates the presence oi a. J matrix. 

Operation 4: Finally, the only term we did not consider (the one with one resolvent and no 
T variable on the last line of ( ID-ll) ) accounts for doing nothing to a dot. 



Substituting eq. fID-ip for each resolvent is equivalent to performing Operations 1 to 4 on all 
dots/links and in all possible ways. After this we have diagrams with free links and missing dots. 
The final step is to join the remaining free links. This is equivalent to contracting the T and T"*^ 
variables in all possible ways. The gluing of free links gives a trivial weight, so this final step is 
just graphical. 

Let us illustrate this procedure on the example treated in subsection 16.21 The following equa- 
tion represents the application of Operations 1 to 4 in all possible ways, with their corresponding 
weights. 







+(1 + 6) 



1 1 

X — a y — i3 

1 1 



1 1 

X + a y + /3 



-O- 



rO : 

A V 



— y"^ — (3"^ ■ ■[[' 







1 1 



1 1 



X — ay — (3 -hD-' x + ay + (3 L-o- 
1 1 



x"^ — a^ y"^ — (3'^ 



(D-2) 



In this equation shaded dots and links represent erased dots and links. Notice that this inter- 
mediate formula can be identified term by term with formula (I6-I2p . Gluing of the free links 
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gives 




x — ay — j3 x + ay + l3 



— oP- -xp- — P"^ 



\ x — ay — p x + ay + pj lq^ 

V^y + p) ( 

0* 



X — a y — p J \ X 



x2 — y"^ — P"^ 

which matches exactly the result (16-130 in section [672] but calculated here using the graphical rules 
we have defined. 

D.2 Computation of the weight for the final diagrams. 

Consider a reduced problem where no erasing of dots is allowed, i.e. only Operation 3 and 4 
are taken into account. In this case only cutting and gluing is allowed and no difference at the 
graphical level appears between the J-antisymmetric and the J-antisymmetric cases. 

To each link in the final diagram G' is attached a weight coming from the different ways we 
obtain it from the original diagram G. That is, when a link in G' is part of G, we can either cut 
that original link and glue it again, or just do nothing. Instead, if the link in G' does not belong 
to G, the only way to obtain it is by gluing cut links. Both contributions will add up to 

+ if the link e G' 

if the link ^ G" 



G G 



G G 



G G 



G G 



1 



1 



1 



Xi-a Vj-P 

1 



Xi~a Vj-fS 
1 ' 1 



1 1 



I 1 1 



1 1 



if the link G G' 

if the link ^G' 

if the link G G' 

if the link ^ G' 

if the link G G' 

if the link ^ G' 



Xi-a Vj+P 

Returning now to the original complete problem, where the erasing of dots is allowed, one 
notices that the weight found in Operations 1 and/or 2 by erasing a pair of dots is the same as the 
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one obtained by forming a cycle with this pair of dots with Operations 3 and 4. Indeed, consider 
a minimal cycle {i.e. a cycle of length 2) in G' and erase from G the dots in this minimal cycle 
using Operation 1 and if possible Operation 2. Operation 1 contributes a factor (1 +6); if the 
minimal cycle in G' is non twisted, we assign it the weight 1, while if it is, we assign it the weight 
Iterating this operation for all minimal cycles in G", one finds that the two procedures, erasing 
pairs of dots or forming minimal cycles with the same dots, produce the same weight. A slightly 
different manipulation is needed when the minimal cycle is present already in G. The outcome 
will be the same. 

Here is an example of the kind of diagrams which have the same weight within the complete 
problem: 

/ 





where A^(G") = M.q is the weight associated with G', for G the original diagram. 

We proved that for the complete problem, the recursion matrix is the same as that of the 
reduced problem. Reassembling everything together and using back the sets of signs s and s', we 
get the recursion matrix 



M 



■G' 
G 



n 

.((x,,s(i)),(x,(,),s'(7r(i))))GG' 



1 



n 

((x„s(i)),(?/.(.),s'(7r(i))))0G" 



n 

.((x.,s(i)),fo,/(,),s'(vr'(i))))eG' 



s{i)xi + as'{'K{i))y^(i) + [3 

1 1 

s{i)xi + as'{Ti{i))y^{i) + (3 

1 1 



1 + 



n 

((x„s(i)),(y^,(,),s'(7r'(i))))^G' 



s{i)xi - a s'{'K'{i))y-^i{i) - (3 

1 1 

\s{i)xi - a s'(7r'(i))?/^/(i) - p 



^This is the origin of the additive ±l's coming with every minimal cycle in the basis of functions eq. (|6-18p . 



52 



Labelling each diagram using the labels in eq. (16-181) . eq. fl6-19p and eq. (16-221) we get 
which completes the proof. 

E Relations between orthogonal/symplectic and unitary 
recursion equations 

In this appendix we relate tetrads uj = {a,T,s,t} introduced in sect. 16.3.11 and permutations 
G &2R, and more precisely to show the bijection between the set of equivalence classes [u] 
and &2R] this leads to an important relation between the recursion matrix Ai and the basis of 
correlation functions F'^ for the orthogonal/symplectic case and the recursion matrix Ai 

-M?;HW,M,«,/?) = li (^^wxw + (E-1) 

\ Xi i/7r(j) y / 



and the basis of correlation functions F 



u 



found in [12] in the unitary case. 

E.l Bijection between <S2R and equivalence clases in GrxQbX x 

Consider pairs of permutations a and r belonging to &r, pairs of sets of R signs s and t belonging 
to As explained in sect. 16.3.11 aoT~^ represents a permutation of R (black) points with signs 
s{i) attached to them, and r o a permutation of R (white) points with signs t{i). To get one 
representative of the equivalence class [{a, r, s, t}], we fix one sign Sj in every cycle of a o r^^ to 
be +1. 



Lemma E.l There is a bijection between &2R (^nd the equivalence classes of Gr x x x Z^. 
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Proof: To construct the bijection we take a permutation vr G &2r- We will relabel the indices 
i = 1, . . . , 2R and call them a{i) = 1, . . . , R, —R, . . . , — 1 (see figure IKTI) . 



i a i<R 

i-{2R+l) iii> R ' 



Define also the auxiliary "rainbow" permutation e{i) = 2R — i for which = id and e{a{i)) = 
«(e(i)) = —a{i). 

12 3 R -R -3 -2 -1 



12 3 



R R+1 2R-2 2R-1 2R 



a 





o o o o o o o o 



12 3 
12 3 



R R+1 2R-2 2R-1 2R 

R -R -3 -2 -1 



a 




Figure 1: Representation of the bijection. The set of lines, irrespective of their type (solid or 
dashed) represents vr G &2r- Solid lines represent a and dashed ones represent r. The signs s, 
resp t, are the signs at the origin, resp.the end, of the solid lines. The arrows, though redundant, 
are meant to help the reader follow the iterations 1-4 above. 



With this new labelling we construct the permutations a and r and the signs s and t as follows 

• 1. Begin with i = 1. 

• 2. Set s(|a(z)|) = sgn(a(z)), then j := 7r(z) and set (T(|a(z)|) = \a{j)\ 

• 3. Set t(|a(j)|) = sgn(a(j)), then k = en~^e{j) and set r~^(|a(j)|) = \a{k)\ 

• 4. If we do not close a cycle set i = k and go back to 2. 
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When we close a cycle of a o r~^, (for example, to close the first cycle we must find again k = 1 
at the end of step 3.), we must open a new cycle. To do so, look at which positive a-type indices 
we have not used yet and choose the smallest one. Set i equal to this value and restart from point 
2. When there is no positive index left, the last cycle is completed and the tetrad {a,T,s,t} is 
constructed. □ 

This procedure is illustrated in figure [ETTI there, the permutation in &2R is tt = (l,2,i? + 
2, 2R, R + 3, -R + 1, 3)(-R)(- ■ ■ ). Following the above rules we determine from it a = (1, 2, R, 3)(- ■ ■ ), 
r = (l,3)(2)(i?)(---), s = (+,-,+,...,-) and t= (+,+,-,...,-). 

Note that the first s sign of every cycle is positive by construction. Note also that conversely, 
constructing tt from the tetrad {a, r, s, t} can be done with the same kind of procedure: this is 
clearly seen in the example. By construction, in this reverse operation, vr depends only on the 
equivalence class [u]. Since the method is deterministic in both directions we have a bijection. 

E.2 Relation between M^^^ and A^^^i?) 
Consider now the set of variables 

{x}2R= {Xi, . . . ,Xr,Xr+i = -Xr, . . . ,X2R = -Xi} , {x}r = {xi, . . . , Xr} 

{y}2R= {yi,...,yR,yR+i = -yR,...,y2R = -yi}, {y}R = {yi, ■ ■ ■ ,yR} 
and consider also the pairs of indices 



where the indices are related through the bijection shown above. Using these definitions it is easy 
to verify that 



TT, it' e &2R 



[{a, a', s, s'}] , [{r, r', t, t'}] e Equivalence Classes of &r x &r x x Z2 



(E-3) 



-jjiR) {T,r',t,t'} 
{t,'t',s,s'] 



■.'}{{x}R,{y}R,a,P) 



X 




5<j'(j)r'(j) Ss(i)t{i) Ss'{a' {i))t' (a'(i)) + 



s{i)xi + a s'{a{i))y„(i) + (3 



1 1 



1 1 



) 




2R, -a, -f3) 
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i.e. the two matrices encountered in the orthogonal/symplectic and unitary cases are in fact the 
same. In this calculation, we have used the bijection to reexpress the Kronecker 5 symbols: 

^a(i),T{i)8s(i),t{i)5s'{a(i)),t'{a(i)) = <57r(Q-l (s(i)i)),7r'{Q-l (s(i)i)) 
Sa'[i)T'[i)Ss(i)t{i)Ss'{a'(i))t'[a'[i)) = <^7r(a-l(-s(j)i)),vr'(a-l(-s(i)i)) 

and to relate the {x}r and {y}R variables with the {x}2r and {y}2R variables according to 
E.3 Relation between F'^ and 

Finally consider the basis of correlation functions in the unitary case for 2R X-type and 2R 
y-type resolvents. In particular consider the components F^^^^{{x}2r, {y}2R, A, B). Call [u] = 
[{(J, r, s, t}] the equivalence class of tetrads corresponding to vr. 

The function F^^^^{{x}2r, {y}2R, A, B) can be constructed from the kind of diagrams shown 
in figure IE. II by following the ireir'^e cycles. These cycles are, by construction, the same ones 
we follow with cr'^, with the only difference that vr is a permutation of 2R elements instead of 
R. Because of this, 7re7r~^e contains two different representatives of the equivalence class [u], and 
since the functions F^^ are independent of the class representative, F'^ contains twice the same 
function F'^ . We can write this as 

F^^,^,{{x}2R,{y}2R,A,B) = {F;^{{x}R,{y}R,A,B))\ (E-4) 

The sign in eq. f l6-22p is, however, not easy to read from vr and e. 
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